quietly{
clear
*Set directory for sampling input and output
cd "/Users/martin/Dropbox/Documents/ProjectMultilateralAid/Data/Rework2021/ReplicationWBER"
use "sampdat.dta"

*first stage
dirifit DR_sec_share_WB_f_11-DR_sec_share_WB_f_74 if DR_sec_share_WB_f_11!=., alphavar(WB_rec_yr idatotal reflow_wbtot libor12_fy ida_WBry refl_WBry lib_WBry) alpha1(WB_secshare_11 WBry_WBsec1 ida_WBry_WBsec1) alpha2(WB_secshare_12 WBry_WBsec2 ida_WBry_WBsec2) alpha3(WB_secshare_13 WBry_WBsec3 ida_WBry_WBsec3) alpha4(WB_secshare_14 WBry_WBsec4 ida_WBry_WBsec4) alpha5(WB_secshare_15 WBry_WBsec5 ida_WBry_WBsec5) alpha6(WB_secshare_16 WBry_WBsec6 ida_WBry_WBsec6) alpha7(WB_secshare_21 WBry_WBsec7 ida_WBry_WBsec7) alpha8(WB_secshare_22 WBry_WBsec8 ida_WBry_WBsec8 ) alpha9(WB_secshare_23 WBry_WBsec9 ida_WBry_WBsec9) alpha10(WB_secshare_24 WBry_WBsec10 ida_WBry_WBsec10) alpha11(WB_secshare_25 WBry_WBsec11 ida_WBry_WBsec11) alpha12(WB_secshare_31 WBry_WBsec12 ida_WBry_WBsec12) alpha13(WB_secshare_32 WBry_WBsec13 ida_WBry_WBsec13) alpha14(WB_secshare_33 WBry_WBsec14 ida_WBry_WBsec14) alpha15(WB_secshare_41 WBry_WBsec15 ida_WBry_WBsec15) alpha16(WB_secshare_42 WBry_WBsec16 ida_WBry_WBsec16) alpha17(WB_secshare_43 WBry_WBsec17 ida_WBry_WBsec17) alpha18(WB_secshare_51 WBry_WBsec18 ida_WBry_WBsec18 ) alpha19(WB_secshare_52 WBry_WBsec19 ida_WBry_WBsec19) alpha20(WB_secshare_53 WBry_WBsec20 ida_WBry_WBsec20) alpha21(WB_secshare_70 WBry_WBsec21 ida_WBry_WBsec21) alpha22(WB_secshare_72 WBry_WBsec22 ida_WBry_WBsec22) alpha23(WB_secshare_73 WBry_WBsec23 ida_WBry_WBsec23) alpha24(WB_secshare_73 WBry_WBsec24 ida_WBry_WBsec24)

predict pxb1, equation(#1) xb
predict pxb2, equation(#2) xb
predict pxb3, equation(#3) xb
predict pxb4, equation(#4) xb
predict pxb5, equation(#5) xb
predict pxb6, equation(#6) xb
predict pxb7, equation(#7) xb
predict pxb8, equation(#8) xb
predict pxb9, equation(#9) xb
predict pxb10, equation(#10) xb
predict pxb11, equation(#11) xb
predict pxb12, equation(#12) xb
predict pxb13, equation(#13) xb
predict pxb14, equation(#14) xb
predict pxb15, equation(#15) xb
predict pxb16, equation(#16) xb
predict pxb17, equation(#17) xb
predict pxb18, equation(#18) xb
predict pxb19, equation(#19) xb
predict pxb20, equation(#20) xb
predict pxb21, equation(#21) xb
predict pxb22, equation(#22) xb
predict pxb23, equation(#23) xb
predict pxb24, equation(#24) xb

*generate predicted values
gen E_WB_share_1=exp(pxb1)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_2=exp(pxb2)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_3=exp(pxb3)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_4=exp(pxb4)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_5=exp(pxb5)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_61=exp(pxb6)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_7=exp(pxb7)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_8=exp(pxb8)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_9=exp(pxb9)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_10=exp(pxb10)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_11=exp(pxb11)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_12=exp(pxb12)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_13=exp(pxb13)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_14=exp(pxb14)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_15=exp(pxb15)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_16=exp(pxb16)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_17=exp(pxb17)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_18=exp(pxb18)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_19=exp(pxb19)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_20=exp(pxb20)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_21=exp(pxb21)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_22=exp(pxb22)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_23=exp(pxb23)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))
gen E_WB_share_24=exp(pxb24)/(exp(pxb1)+exp(pxb2)+exp(pxb3)+exp(pxb4)+exp(pxb5)+exp(pxb6)+exp(pxb7)+exp(pxb8)+exp(pxb9)+exp(pxb10)+exp(pxb11)+exp(pxb12)+exp(pxb13)+exp(pxb14)+exp(pxb15)+exp(pxb16)+exp(pxb17)+exp(pxb18)+exp(pxb19)+exp(pxb20)+exp(pxb21)+exp(pxb22)+exp(pxb23)+exp(pxb24))

*second stage
dirifit f_sec_share_bilat_11 -f_sec_share_bilat_74  ,  alphavar(ln_pop ln_gdp l_committed_bilat_yr l_WB_rec_yr) alpha1(sec_share_bilat_11 E_WB_share_1) alpha2(sec_share_bilat_12 E_WB_share_2) alpha3(sec_share_bilat_13 E_WB_share_3) alpha4(sec_share_bilat_14 E_WB_share_4) alpha5(sec_share_bilat_15 E_WB_share_5) alpha6(sec_share_bilat_16 E_WB_share_6) alpha7(sec_share_bilat_21 E_WB_share_7 ) alpha8(sec_share_bilat_22 E_WB_share_8) alpha9(sec_share_bilat_23 E_WB_share_9) alpha10(sec_share_bilat_24 E_WB_share_10) alpha11(sec_share_bilat_25 E_WB_share_11) alpha12(sec_share_bilat_31 E_WB_share_12) alpha13(sec_share_bilat_32 E_WB_share_13) alpha14(sec_share_bilat_33 E_WB_share_14) alpha15(sec_share_bilat_41 E_WB_share_15) alpha16(sec_share_bilat_42 E_WB_share_16 ) alpha17(sec_share_bilat_43 E_WB_share_17) alpha18(sec_share_bilat_51 E_WB_share_18) alpha19(sec_share_bilat_52 E_WB_share_19) alpha20(sec_share_bilat_53 E_WB_share_20) alpha21(sec_share_bilat_70 E_WB_share_21) alpha22(sec_share_bilat_72 E_WB_share_22) alpha23(sec_share_bilat_73 E_WB_share_23) alpha24(sec_share_bilat_74 E_WB_share_24) robust cluster(r_ccode) 
*NOTE: can't use time-smoothened IVs bc time dimension gets lost in resampling. 

*save coefficients, turn into variables
matrix b=e(b)
svmat b, names(b)

keep b*
drop bd_WDI
duplicates drop
drop in 2

save "sampout.dta", replace
}

